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ABSTRACT 

We consider a supersymmetric model with a single matter supermultiplet in a five-dimensional 
space-time with orbifold compactification along the fifth dimension. The boundary conditions 
on the two orbifold planes are chosen in such a way that supersymmetry remains unbroken on 
the boundaries. We calculate the vacuum energy-momentum tensor in a configuration in which 
the boundary branes are moving with constant velocity. The results show that the contribution 
from fermions cancels that of bosons only in the static limit, but in general a velocity-dependent 
Casimir energy arises between the branes. We relate this effect to the particle production due 
to the branes motion and finally we discuss some cosmological consequences. 
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1 Introduction 



The presence of non-trivial boundary conditions in quantum field theory leads to interesting 
phenomena such as the shift in the zero-point energy with respect to the unbounded space 
configuration [1]. In the case of electromagnetic fields between perfectly conducting parallel 
plates, such shift has been observed as an atractive force acting on the plates (see [2] for 
a recent measurement). Since the original work of Casimir, this effect has been studied in 
different contexts, thus for instance, it plays a fundamental role in hadron physics, in the 
so-called bag model in which quarks and gluons are confined inside a spherical shell and the 
corresponding Casimir energy contributes to the total hadron mass. Also in higher-dimensional 
theories, the stability of the compactified extra dimensions depends crucially on this effect. 
Other applications can be found in condensed matter physics and atomic or molecular physics 
(for a review see [3]). 

More recently, the increasing interest in the construction of cosmological models based on 
string or M theory has focused on the so-called brane-world scenario, in which our universe is 
understood as a 3-brane living in a higher dimensional space-time. A particularly interesting 
model appears in the low-energy regime of M-theory described by eleven dimensional super- 
gravity [4]. After compactification of six dimensions in a Calabi-Yau manifold, we are left with 
a five-dimensional model in which matter fields live on the two boundary branes arising after 
S* 1 /^-orbifolding the fifth dimension. In this context, the Casimir energy arising between the 
two static boundaries has been computed in [5, 6], in the first of these two works, the back- 
reaction on the geometry was taken into account. The same problem has been considered in 
five-dimensional anti-deSitter space in [7]. The recently proposed ekpyrotic (cyclic) model of 
the universe [8] is also based on this framework in which the motion and collision of two such 
branes is responsible for the Big-Bang of the standard cosmology. 

In this paper we are interested in studying the possible effects of the Casimir energy in an 
scenario like the one mentioned before in which two branes are moving towards each other. The 
complete analysis of the problem is in general too involved to obtain explicit analytic results 
and, for that reason, we will consider a simplified model in which the two branes are moving 
with constant relative velocity and they are perfectly flat, ignoring possible gravitational effects. 
In any realistic model of a brane collision process it will be necessary to consider the acceleration 
and the brane curvature [9], but the present analysis would be the first (velocity-dependent) 
correction to the flat static case. We will also consider a simple global supersymmetric model in 
five dimensions with a single matter hypermultiplet instead of working with the full supergravity 
lagrangian in eleven dimensions. Despite these approximations the results will shed light on the 
velocity dependence of the Casimir energy and therefore on the stability of the branes system. 
We will be interested also in the role of supersymmetry in the cancellation of the Casimir energy 
in the moving case. For that purpose we will need to extend previous studies of the dynamical 
Casimir effect [10] to fermionic fields. In fact we will show that only in the static limit the 
contribution from fermions cancels that of bosons (as expected from supersymmetry), but in 
general a velocity dependent Casimir energy arises between the branes. 

The paper is organized as follows, in Section 2, we introduce the five-dimensional super- 
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symmetric model and set the boundary conditions for the different fields. In Section 3 we 
compute the scalar contribution to the Casimir stress tensor using the image method for Green 
functions for arbitrary value of the brane velocity. In Section 4, we perform a similar calcu- 
lation but with the Bogolyubov transformations technique in the low-velocity regime. Section 
5 is devoted to the fermionic contribution, and finally Section 6 contains a discussion on the 
possible cosmological consequences, paying special attention to the ekpyrotic (cyclic) model of 
the universe. 



2 Supersymmetric model Lagrangian in 5-D 

Let us consider the simplest globally supersymmetric model in 5-D. A five-dimensional off-shell 
hypermultiplet consists of two complex scalar fields, which we denote by 0i and <f) 2 , one four- 
component Dirac spinor ip and two auxiliary complex scalars fields F and G [11, 6, 12]. The 
smallest representation of the Clifford algebra in five dimensions is four dimensional and is given 
by: 7 A = {7 , 7 1 , 7 2 , 7 3 , 7 4 } where 7 4 = ry 5 with the usual definition 7 s = i'j ^ 1 ^ 2 "/ 3 . With 
this definition we have {^ A ^ B } = 2i] AB where r\ AB = diag(+, — , — , — , — ). As a consequence 
the chirallity operator 7°7 1 7 2 7 3 7 4 = 1 is trivial and it is not possible to define chiral fermions 
in five dimensions. 

The simplest supersymmetric action corresponds to massless non-interacting fields and is 
given by: 

S = Jd 5 x (r) AB dA<t>ld B cf>i + V AB dA<t>*2d B 4> 2 + i^ A d A ^ + \F\ 2 + \G\ 2 ) (1) 

As commented in the introduction, we will assume that the fifth dimension is compactified 
on the orbifold S 1 /Z 2 , that is, the fifth dimension is just an interval [0, L], with two parallel 
boundary planes at x 4 = 0, L. In fact, in order to perform the orbifold projection, we need to 
assign Z 2 parities rj = ±1 to the different fields in such a way that $(x M ,x 4 ) = r]&(x^, —x 4 ) 
where $ is any of the hypermultiplet components. A consistent assignment which respects 
N = 1 supersymmetry on the boundaries is given in Table 1 (see [13, 6]), where we have used 
the definitions -0L = Pi$ = ^f-^P an d ipR = Pr^P = 

For an S* 1 compactification with radius R, the Fourier expansions for the fields with different 
Z 2 parities become: 

00 

$+(^,X 4 ) = J2 C0S ( nX V R )®n( xfl ) 
n=0 

00 

x 4 ) = sin(nx 4 / J R)$ n (^) (2) 

n=l 

We see that even-parity fields satisfy Neumann conditions on the boundary planes, whereas 
odd-parity fields satisfy Dirichlet conditions i.e. for L = irR we have: 

d^ + (x»,x 4 ) \ x 4 =0jL = 0, 

$_(^,x 4 )|, 4=0ii = 0. (3) 
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1] = 1 
(Neumann) 


i] = — 1 
(Dirichlet) 


01 


02 






F 





Table 1: Parity assignments for the model in (1) 



If supersymmetry remains unbroken we expect that the contribution from fermions cancels 
that of bosons in the vacuum energy-momentum tensor. As we will see this is the case when 
the boundaries are static, however here we will consider a more general situation in which the 
boundaries are moving. 

3 Scalar contribution 

Let us write the classical energy-momentum tensor for massless real scalars in a five-dimensional 
Minkowski background: 

T S AB = (1 - 2O0>,b + (2f - 9ab9 CD <P,c<I),d 
2 

- 2£0,AB0 + -77AS0D0 (4) 


where for generality we have considered arbitrary non-minimal couplings in the scalar lagrangian 
(,R(fi 2 , which vanish in flat space but however have a no n- vanishing contribution to the energy- 
momentum tensor. The particular case £ = 3/16 corresponds to the conformal coupling in five 
dimensions 

In order to calculate the vacuum expectation value of the energy-momentum tensor we will 
use the Green's function method. With that purpose we introduce the so called Hadamard 
function: 

D^\x,y) = (O|0(:r)0(y) + 0(y)0Or)|O> (5) 

where |0) denotes the usual Minkowski vacuum state built out of the plane- wave solutions of 
the free equations of motion written in cartesian coordinates. It is possible to write the vacuum 
expectation value of Tab in cartesian coordinates as: 

(0|Tf s |0) = , lim \ ((1 - 20ajag + - I) VABV CD d x cd y D - 2idfdi) D^(x',y') (6) 
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where we have used O x 'D^(x' : y') = 0. In the simplest case, i.e., unbounded Minkowski space, 
the Hadamard function can be obtained from the mode expansion of the quantum fields and it 
is given by: 



1 e(a(x,y)) 



±iy)l ( 7) 

4vr 2 a(x,yf/ 2 1 ' 

where a(x,y) = —tiab{.x a — y A ){x B — y B ), 9(a(x,y)) is the step function and d 4 p is the inte- 
gration measure for the spatial components of the five- momentum p A - 

We will first consider the case in which the two branes at x 4 = 0, L are fixed. The Hadamard 
functions satisfying Dirchlet and Neumann conditions can be calculated by the image method 
(see for example [13]). According to the principle of mirror reflection such functions can be 
obtained as a sum over an infinite number of the free functions in (7) evaluated at the image 
points x' n . They can be written as: 

oo 

D { n) N M= £ (DW(x,x> n+ ) T DW(x,x> n _)) (8) 

n=— oo 

where the minus sign corresponds to the Dirichlet function and the plus sign to the Neumann 
one. The image points are given by: 



x 



n± 



(t',x'\x' 2 ,x ,3 ,±x >4 -2nL) (9) 



These functions satisfy the corresponding boundary conditions given in (3). Since in our model, 
one complex scalar field satisfies Dirichlet conditions and the other one Neumann conditions, the 
total scalar energy-momentum tensor will be given by the sum of the two kinds of contributions, 
i.e.: 



<0|T? B |0>= hm 2 ((1 - 209^ - l&i&i 

x,y'—>x \ 

+ - i) VABV CD d x cd£) £ D( V, y' n+ ) (10) 

The expression in (10) is divergent in the coincidence limit, the divergent contribution coming 
from the n = term which is nothing but the unbounded space vacuum energy-momentum 
tensor. Thus we can define the renormalized vacuum expectation value substracting this free 
space term: 

(0\T AB \0) ren = (0\T AB \0) - (0\T AB \0) n=0 (11) 

The final result does not depend on the non-minimal parameter £ and is given by the well-known 
expression: 

(0|Ti B |0)^ c = -|^j|o -1 | (12) 
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Let us consider next the case in which the 3-brane at x A = is fixed whereas the other one 
is moving with constant velocity v in the fifth dimension. In order to obtain the image points it 
is easier to work in a coordinate system in which the branes are fixed, but the geometry is con- 
tracting (expanding) along the fifth dimension [14]. Thus, let us consider the new coordinates 
3 , x) defined as: 

x° = r cosh(x), x 4 = rsinh(x) (13) 

In these ccordinates the branes positions are given simply by x — and x — Xo- Notice that 
curves with constant fifth coordinate x — Xo describe brane motions with a constant velocity 
v = tanh(xo)- The corresponding metric tensor reads: 

rfs 2 = dr 2 _ ( dx iy _ ((fa 2 )2 _ ^3)2 _ T 2 dx 2 (14) 

which is nothing but the Milne metric in five dimensions. This metric can describe the portions 
of Minkowski space-time with either x° > or x° < 0. We will use indices A,B,C,... = 
0, 1, 2, 3, 4 for Minkowski coordinates and M, N, P, ... = 0, 1, 2, 3, 4 for the Milne coordinates. 
For this metric the only non- vanishing Christoffel symbols are Y T XX = a 2 T, T* x = T* T — K In 
the following, we will use the hat to denote objects written in Milne coordinates. The invariant 
interval a(x,x') defined before can be written in the new coordinates as: 

a(x, x) = -t 2 - t' 2 + 2tt' cosh(x - x) + (x 1 - x' 1 ) 2 

+ ( x 2 _ x y + ( x 3 _ x y (15) 

The Dirichlet and Neumann functions are still given by (8), but with the image points written 
in the new coordinates, i.e.: 

%n± = (T',x n ,x' 2 ,x' 3 ,±x - 2nxo) (16) 
The boundary conditions satisfied by the Green functions now read: 

d d\x,x')\x,x'=o,x = 0, 
d x D$(x,x')\ x , x , =0tX0 = (17) 

Using Eqs. (6), (11) and changing to Milne coordinates, we obtain the renormalized expres- 
sions for the energy-momentum tensor. In these coordinates the tensors are diagonal and do 
not depend on the spatial coordinates, however unlike the static case, they do depend on the £ 
parameter. Thus for minimal coupling £ = we get: 



(0|^MJv|0)ren' — or* 2UI5 X! 



(I \ 



327r 2 |r|5^i | sinh 5 (n X o)| 



Q _ 3cosh(2nxo)-l , j 

2 



(18) 



V 4t z J 



whose trace is non-vanishing. Notice that when the branes are fixed the result (12) is traceless 
regardless of the value of the £ parameter. For the conformal coupling £ = 3/16, we get: 
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which is traceless as expected. These results are valid for arbitrary values of the brane velocity. 
However, in order to extract an explicit velocity contribution to the Casimir energy-momentum 
tensor, we will compare them in the non-relativisitic limit t; C 1 with the static result in (12). 
In Milne coordinates, the proper distance between the two branes is given by d = t X o, thus, in 
such limit, we get for minimal coupling: 



(0|T^|0)— * 



3C(5) 

327T 2 rf 5 



/ 1 












4r 2 



\ 



5((3)v 2 

647T 2 rf 5 



1 





13- 



fl 
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\ 


4r 2 J 



(20) 



and for conformal coupling: 



\T S \r\\ mov ~ 

\ 1 MN\ v /ren ~ 



327T 2 rf 5 



1 











-1 











4r 2 



(21) 



The velocity contribution to the Casimir energy is quadratic and positive in the non-relativisitic 
limit. In the following we will show that such contribution is nothing but the energy density 
in the form of particles created from the vacuum by the moving boundaries. 



4 Particle production 

Let us consider a quantized real scalar field: 



^) = EEte n + i M L) ( 22 ) 

% n 

The modes u^ n obeying the Klein-Gordon equation in Milne coordinates are given by: 

u^Jx) = Nf kn (r)e ikS cos(v X ), Neumann 
u£ n (x)=Nf kn (T)e ins sin(v X ), Dirichlet (23) 

where the functions fk n (j) satisfy the equation: 

(dodo + k 2 + ^ + ±A>) fUr) = (24) 
with k 2 = k 2 , v = nir/xo an d n a positive integer number. The solutions of equation (24) can 

(1 2) 

be written in terms of Hankel functions fkn = H\v (^ r ) where the 1(2) functions behave as 
positive frequency modes in the r — > +(— )oo limit. The solutions are normalized with respect 
to the scalar product: 

(0i, <h) = ~i [ d*x r d X r{<P\do<p2 - 40i0 2 ) (25) 
Jv Jo 
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where we are working in a box of finite volume V in the ordinary three-space. We will take the 
continuum limit at the end of the calculations. From the expression above, we can obtain the 
normalization constant iV = \J^^ e ^ vlT ^ 2 - 

Because of the moving boundaries, a given solution which behaves as positive frequency at 
a given time will become a linear superposition of positive and negative frequency modes at 
any later time. This implies that the initial vacuum state will contain particles as seen by an 
observer at a later time. Thus consider the solutions at a given time for the static problem 
with a fixed brane separation L. The positive frequency modes satisfying the different boundary 
conditions are given by: 



v 



f n = Ne -Mkx cos ( gz ) ? Neumann 



vf n = Ne- iuJt+ikS sm(qz), Dirichlet (26) 

where q = nn/L, uj 2 = k 2 + q 2 and iV = (LVuj)^ 1 / 2 . These modes will define a different Fock 
space which can be built out of the corresponding creation and annhilation operators which we 
will denote by 6^ n , frt^. Thus, for either the real or the imaginary part of our scalar fields we 
have: 

^W = EEte n + te (27) 

% n 

The two sets of operators are related in general by the following Bogolyubov transformation 
[15]: 

b kn = E E { a kk>nn< a k>n> + %W fl L) ( 28 ) 
k' n ' 

The Bogolyubov coefficients relating both solutions are given by: nn i = { u k n ^ v p n i)- fol- 
lowing [14] , we are interested in the calculation of the total number of particles created by the 
brane motion in the non-relativisitic limit Xo ^ 1> i- e - f° r large values of v. In that limit, 
we can use the large-order expansion of Hankel functions in the calculation of the Bogolyubov 
coefficients. In the present case it is enough with the first term in the expansion which is given 
by [16]: 



H - ik ^ - f.j^hr exp (t + 1(1/2 + fcV)V2 " ivaicsinh (£) - *i) (29) 

This expression is valid for large r, i.e. large brane separations and satisfies the equation 
of motion up to terms of 0{y~ 2 \ In the non-relativistic limit, we have t = x° = t + 0(x 2 ), 

z — r X + C(x 3 ) an d v = Xo + CKXo)) so th & t we can write L = txo + C(Xo) an d q = v /t + O(xo) 
so that lu 2 = k 2 + u 2 /t 2 + 0(xl)- Therefore for the Bogolyubov coefficients relating Neumann 
solutions we get: 



fe„n' " Jy d3x l dZC0S 



L . (n'irz\ ( nnx\ e i{k ' k)£ k 2 r 

L J 00 V XO J 2 X0 (z/2 + A ;2 r 2)3/2 e 
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(30) 



where the phase is given by: 

$(z/, k, t) = 2ur - z/arcsinh (j^J + ^ (31) 

and the same results is obtained for Dirichlet solutions, whereas mixed Dirichlet-Neumann 
coefficients vanish. Defining /%d / = fl^kk'^nri, the total energy-density of the particles 
created for the scalar matter content in our model with two complex scalar fields is: 

ps = l^J (2^H^' * z n ? 7 (32) 

where we have already taken the continuum limit V — > oo. This is a divergent expression which 
can be renormalized using the zeta-function method [17]. By means of a change of variables, 
we can write the above expression as: 

ps ^^p%L dx w 2 +» 2 r /2 = T6^- B {l> - 1 ) S (^) (33) 

where B(p,q) = T(p)T(q)/T(p + q) is the beta function and J2n n ~ s = C( s ) is the Riemann zeta 
function. We can cancel the divergent contributions from the integral and the series by using 
the analytical continuation of the relation: ((s) = T 71-5 X ^ 2 C(1 ~~ S )/T (§), we obtain the 
following renormalized result: 

Ps * a A 2 5 5 5 C(3)^ 2 (34) 

which agrees with the first correction to the static Casimir energy in (20). We see that this 
velocity correction can be understood as the contribution from the particles produced by the 
non-adiabatic motion of the branes [14]. 



5 Fermionic contribution 

Although the fermionic Casimir effect has been considered previously in different contexts, only 
a few results in the case with moving boundaries can be found in the literature [18]. In this 
section we will consider the problem in our five dimensional model. 

The fermionic energy- momentum tensor for massless fermions in five dimensions can be 
written in Minkowski coordinates as: 

t ab = 7 (^7a<9 b V + ^IbOa^ - Oa^Ib^ - d B ^lA^) (35) 
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Following similar steps to the scalar case, we will compute the corresponding vacuum expecta- 
tion value, in the case in which the two branes are fixed, by means of the image method. The 
fermionic Hadamard function can be written as: 

= {o\My)Mx)-Mx)My)\<>) (36) 

where a, b are spinor indices. From (35), the vacuum expectation value of T^ B can be written 
in Minkowski space-time in terms of the Hadamard function as: 

(0|TL|0) = , lim 1 {{ lA ) ba di + ( lB ) ba di) V$(x>, y') (37) 



In the simplest case with no boundary conditions and in Minkowski coordinates, the Hadamard 
function is given in terms of the scalar function (7) as V^ b \x,y) = — i^ b dj^D^(x, y). In the 
case with fixed boundaries, the different fermionic components satisfy (3) and the Hadamard 
function obeying these boundary conditions can be written as: 

VU(x,y) = P L V%(x,y)P R + P R V%(x,y)P L 

+ P L V { S N {^y)P L + P R V%{x,y)P R (38) 

With these relations, we obtain from (37) the following expression for the energy-momentum 
tensor between the fixed branes: 

(0|Tj» = , lim \tr(diV%{x\y')P RlA P L + diV%{x\y')P LlA P R 

x'—*x, y'—>x 4 v 

+ d x B V%(x', y')P LlA P L + ViV%(x\ y')P RlA P R ) + (A^B) 

(39) 



In order to calculate the explicit expression of the different components of the Hadamard 
function in (38), we need to know the mode expansion of the Dirac field: 

m = + ( 4 °) 

kna 

where we have labelled the states by a = ±1/2, according to the value of the spin projec- 
tion along the ordinary three-momentum on the brane. The positive and negative frequency 
solutions of the Dirac equation in Minkowski coordinates are given by: 

= ,-„}. . ( IS? ) ( 41 ) 



^2LVu{u ± q) V ^s(qz)U a 



where the two-component spinors with the appropriate helicities are given by: Uy 2 = (1, 0) 
and UT.\/2 = (0) 1) an d the component of the momentum in the fifth dimension q is defined as 
in the scalar case with fixed boundaries. The modes are normalized with respect to the scalar 
product: 

{i>\M = J d^Xy/g^-riQ^ (42) 
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where go denotes the determinant of the induced metric on the hypersurfaces of constant time. 
In the calculation of the energy-momentum tensor in (39) only derivatives of the Green functions 
are relevant, thus, making use of the the equations of motion, ^ x V^\x, y) — and acting with 
the projectors we get: 

1*dZP R V$> D (x, y)P R = -rd*P L V%(x, y)P R (43) 

and a similar expression is obtained for the NN function. Tracelessness and conservation 
determine a diagonal form of the energy-momentum tensor in the Minkowski background. 
Accordingly the knowledge of the DD and NN functions together with Eq. (43) is sufficient 
to compute the complete tensor. Inserting the explicit mode expansion (40) into (38), it can 
be shown that such functions are given by: 

P L V%(x, y)P R = -iP L ^D^\x, y)P R (44) 

where Dp\x, y) is the scalar function given in (8) and a similar expression is obtained for the 
NN function. From (39) we finally get: 

(0|TL|0)f e f c = 2 , lim ± didiD^(x',y' n+ )5 AB = ^- 5 [ -1 ) (45) 

which agrees with the static scalar case (12) but with opposite sign. Accordingly, and as 
expected from supersymmetry, for static branes the total energy-momentum vanishes. 

The case in which one of the branes is moving cannot be computed easily with the image 
method, so that in the following we will use the Bogolyubov transformations technique. Let us 
then consider again Milne coordinates. We introduce the following expression for the vielbein 
e A M and the inverse vielbein corresponding to the Milne metric in (14): = = 5^ u 
with /i — 0, ..,3, and e\ = r, e\ = 1/r. The curved gamma matrices are defined as usual by 
7 M = e^" 4 . The Dirac equation in these coordinates reads: 

iflu = + l 7 4 d 4 + ±^ u = o (46) 




We look for solutions in the form up = i&i>f with: 

kna kna 

where ± index denotes positive or negative frequency and again the fifth component of the 
momentum v is defined as in the scalar case with moving boundaries. With these definitions 
it can be seen that the spinor obeys the correct boundary conditions, i.e.: 

d x P R ulJ x=0M = (48) 
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The function f^Jj) satisfies the equation: 

/!» + (* 2 + % ~ h + %) 4W + ;/f>) = o (49) 

The corresponding solutions are given by Hankel functions f^Jj) = H^_y 2 (kr) and f^ n ( r ) = 

Hiv-i/2\k T )- m ^ ne non-relativistic limit i/ 1 the asymptotic expansion of the positive 
frequency Hankel function reads [16]: 



if« (jfe^ ~ > - 1 1 w \ 

n i"-W KT) ~ V 7T (FT 2 + 1/2)V4 ^ 4(^2 + ^2)^1 

x exp (i(£; 2 T 2 + i/ 2 ) 1 / 2 + ^ - (ii/ - 1) axcsinh (50) 

This expansion satisfies Eq.(46) up to 0{y~ 2 \ The normalization constant can be obtained 
from the scalar product (42) and is given by: iV = ig-W 2 ^ ^^ . Following similar steps to the 
scalar case, we calculate the number of fermions produced by means of the Bogolyubov coeffcient 
%',m>«< = ("jL'^n'J- After a len S th y calculation we obtain %, ;nn>Q , = PknJkk' S nn'S aa ' 
with /3g na = to the lowest order i.e. O(v ) and 



4fcz/ (jfe 2 + £) ; 



x exp ^2iwr — if arcsinh ^^~^) + (51) 

including the first order ((9(z/ -1 )) correction. 

The total energy density to this order is then given by: 

9 00 r rl 3 k 

p F = r E E/tHs^I 2 ( 52 ) 



L a=±l/2n=0 J ( 27r ) 

where the factor of two comes from the fact that antifermions are produced in the same amount 
as fermions. Using again zeta-function regularization to perform the integrals we obtain: 

+ Kl4) +B (l°))S" 2 = -64^3« 3 ) < M > 



Finally combining this result with the static contribution in (12) and taking into account 
that for massless fermions the theory is conformally invariant, the total result for the renor- 
malized energy-momentum tensor up to second order in the velocity is: 

1 



(0&0>— ^7^(C(5)-^; 2 ) ! -1 j (54) 



32vr 2 d 5 



4r 



2 
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We can compute the total result including minimal scalar and fermion contributions, it is 
given by: 

C(3)v* l l 
(0|T^|0)— = (0|T^|0)— +(0|T^|0)— ^-tH^ "T 1 i (55) 



16*** Q g 4r2 



For conformal coupling we get: 

(0|T^|0) 



tot , m „ 17C(3)^ 5 



1287T 2 rf 5 









- 1 








4r 2 



(56) 



We see that the result is non-vanishing for o / 0, i.e. a Casimir effect arises between 
the branes even though we are using the same kind of boundary conditions which preserved 
supersymmetry in the static case. 



6 Conclusions and discussion 

In this work we have considered the dynamical Casimir effect in a five dimensional global su- 
persymmetric model between two moving boundary branes. Taking the same kind of boundary 
conditions as in the static case, the results with moving boundaries are however completely 
different. In fact a non- vanishing negative Casimir energy is generated both for minimal and 
conformally coupled scalars. This induces a new type of velocity dependent potential between 
the branes and the breaking of supersymmetry due to the branes motion. 

The presence of this additional energy density could have interesting cosmological effects. 
Thus for instance, in the context of the ekpyrotic (cyclic) model [8], the universe contracts to 
a singularity and then reexpands. However it was shown in [19] that ultra-relativistic particles 
produced near the singularity could disrupt the cyclic evolution provided their energy density 
reaches the Planckian scale. We have shown that the velocity contribution to the Casimir effect 
can be interpreted as particle production in the Minkowski vacuum and in fact the energy 
density is expected to be very high in the limit of small brane separation. The stability of this 
kind of orbifold singularities has been also studied in [20]. According to their conclusions, the 
introduction of a single particle between the branes causes the collapse of the universe into a 
curvature singularity. 

Concerning the problem of particle production in a Milne metric, in a recent paper [21], it 
has been shown that by matching modes through the singularity for a free theory, it is possible 
to find an appropriate vacuum state in which particles are not produced. However, when 
including time-dependent interactions, particles are generically produced. We understand that 
this should be the case in the presence of the moving boundaries, implying that the energy- 
momentum tensor calculated in this work is physical and should backreact on the space-time 
geometry. (See also [22] for a discussion about the backreaction problem.) 
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